Chiral Lagrangian at finite temperature from the Polyakov-Chiral Quark Model 
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We analyze the consequences of the inclusion of the gluonic Polyakov loop in chiral quark models 
at finite temperature. Specifically, the low-energy effective chiral Lagrangian from two such quark 
models is computed. The tree level vacuum energy density, quark condensate, pion decay constant 
and Gasser-Leutwyler coefficients are found to acquire a temperature dependence. This dependence 
is, however, exponentially small for temperatures below the mass gap in the full unquenched calcu- 
lation. The introduction of the Polyakov loop and its quantum fluctuations is essential to achieve 
this result and also the correct large N c counting for the thermal corrections. We find that new 
coefficients are introduced at 0(p 4 ) to account for the Lorentz breaking at finite temperature. As 
a byproduct, we obtain the effective Lagrangian which describes the coupling of the Polyakov loop 
to the Goldstone bosons. 
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I. INTRODUCTION 

At zero temperature, confinement and spontaneous 
breaking of chiral symmetry emerge as distinct features of 
QCD. This explains the absence of free quarks and gluons 
as well as the observed mass gap in the hadron spectrum 
between the pseudoscalar mesons and the rest of par- 
ticles and resonances. At a given critical temperature, 
lattice simulations predict a deconfinement phase tran- 
sition where chiral symmetry is simultaneously restored 
(for a review see e.g. Ref. JjJ and references therein). 
This remarkable coincidence between the deconfinement 
and chiral phase transitions remains so far unexplained 
from the theoretical side. Nevertheless, there are two 
extreme limits where each of the phase transitions can 
be characterized by an order parameter. For extremely 
light quarks, the quark condensate is used as an order 
parameter for the chiral phase transition where one goes 
from a non-vanishing to a vanishing value across the 
phase transition. In the opposite limit of infinitely heavy 
quarks, the deconfinement phase transition can be char- 
acterized by a breaking of the center symmetry of the 
gauge group and the order parameter is the Polyakov loop 
which evolves from a vanishing value to unity above the 
critical temperature. The real situation for light quarks 
is in between but one still observes sudden changes both 
in the chiral condensate as well as in the Polyakov loop. 
In the present paper we want to address in a quanti- 
tative manner a rather remarkable feature that arises at 
low temperatures in quark models when the chiral flavour 
symmetry and colour gauge center symmetry are jointly 
considered. 
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Besides the existence of a mass gap, a further out- 
standing consequence of spontaneous chiral symmetry 
breaking is that the would-be Goldstone bosons interact 
weakly at low energies and Effective Field Theory Meth- 
ods such as Chiral Perturbation Theory (ChPT) 0, Q 
(for a review see e.g. Q and references therein), can 
successfully be applied in terms of unknown low energy 
constants (LEC's) which cannot be explained on the ba- 
sis of the symmetry alone. This implies neglecting the 
explicit effects of states about the mass gap which in the 
case of two flavours might be identified with the scalar or 
vector meson masses, and limits the maximum energy at 
which standard ChPT may confidently be applied. Ac- 
tually, the bulk of the values of the LEC's can be satu- 
rated by the low energy contribution stemming from the 
exchange of resonances located in the mass gap region. 
At finite, but low, temperatures the physics of QCD is 
believed to consist of a gas of heated hadrons and one 
still expects the dominant role to be played by the pseu- 
doscalar mesonic thermal excitations 0, B ; effects of 
resonances are exponentially suppressed by a Boltzmann 
factor in the mass gap e~ mp l T . Obviously, the applica- 
bility of such an approach requires that there still be a 
mass gap and that confinement still holds. This enti- 
tles, in particular, to consider the LEC's of the Chiral 
Lagrangian as temperature independent couplings and 
finite (low) temperature model independent predictions 
are deduced 

II EES El From this point of view 
finite temperature ChPT provides a strong theoretical 
constraint on the QCD physics well below the phase tran- 
sition. Extrapolations based on ChPT suggest a melting 
of the condensate. However, it is unclear whether this 
vanishing meets the very requirement of a phase tran- 
sition regarding the quark condensate. Moreover, such 
a purely hadronic based description can never account, 
by construction, for the deconfinement phase transition 
expected from lattice simulations 0. 

As already mentioned, the deconfinement phase tran- 
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sition can be characterized for infinitely heavy quarks by 
a breaking of the center symmetry of the colour gauge 
group. The interplay between this center symmetry and 
chiral symmetry requires explicit consideration of quark 
degrees of freedom and can quantitatively be assessed in 
chiral quark models where the spontaneous chiral sym- 
metry breaking is implemented (see e.g. Ref. [l4j for a 
review and references therein) and the mass gap can be 
identified as twice the constituent quark mass M. Chiral 
quark model Lagrangians are invariant under the flavour 
chiral group SUn(Nf) <g> SU^Nf) and have often been 
used to provide some semiquantitative understanding of 
hadronic features in the low energy domain. In addition, 
at zero and finite temperature the standard chiral quark 
models are invariant under global S\J(N C ) transforma- 
tions. At some critical temperature Chiral quark models 
predict already at the one loop level a chiral phase tran- 
sition 0, 0| at realistic temperatures. Pion corrections 
where first considered in Ref. [13 (for a review see e.g. 
Ref. 0). This has been traditionally considered a big 
phenomenological success for these models, but they suf- 
fer from the unphysical contribution of states which are 
not colour singlets, so that even at lowest temperatures 
the hot environment corresponds to a plasma of multi- 
quarks, characterized by Boltzmann factors ~ e - nM / T 
and not colour neutral hadronic states [Tsj. This un- 
physical feature can be avoided by noting the relevance 
of large and local SU(iV c ) gauge invariance at finite tem- 
perature. Basically, this corresponds to include non- 
perturbative finite temperature gluons. Based on previ- 
ous works I2^ . l2lll22| (see also we have discussed 
at length [ljj how large gauge invariance can efficiently 
be implemented by the coupling of the Polyakov Loop to 
the standard chiral quark models hence providing a cool- 
ing mechanism to the chiral-deconfincmcnt phase transi- 
tion and the melting of the chiral condensate is shifted 
up from 150 MeV to about 250 MeV in more quantitative 
agreement with QCD lattice simulations. The coupling 
of QCD distinctive order parameters at finite tempera- 
ture to hadronic properties has been the sub ject of much 
attention over the recent past [2(J, [22, 12a, |2y, [23, l2q| . 
Effective actions for the Polyakov loop as a confinement - 
deconfincment order parameter have been proposed be- 
cause of their relevance in describing the pha se transition 
from above the critical temperature [23, l30l l3ll l33 | . 

In the present paper we want to analyze in more detail 
these Polyakov-Chiral Quark Models and, more specif- 
ically, how the Polyakov loop couples to the lightest 
Goldstone bosons at low energies, and what is the net 
effect at low temperatures in the effective Chiral La- 
grangian. Amazingly, the inclusion of Polyakov loops, a 
colour source, at a quantum level reproduces the expec- 
tations of low temperature ChPT based on the existence 
of a mass gap, namely the tree level coefficients of the chi- 
ral Lagrangian are exponentially suppressed in the mass 
gap. While this supports the ChPT-expected dominance 
of pion fluctuations at low temperatures, the full model 
unlike ChPT still predicts a rapid decrease of the mass 



gap at higher temperatures as well as a sudden rise of 
the Polyakov loop at the same temperature. 

The paper is organized as follows. In Sect. [H] we de- 
scribe shortly the Polyakov-Chiral quark models, and 
more specifically the simplest version of them. In Sect. lIIII 
we set up the calculational framework of the low energy 
chiral Lagrangian at finite temperature and display al- 
ready the general structure of the main result. In Sect. lIVl 
we use the technique of heat kernel expansion at finite 
temperature with Polyakov loops and the associated large 
gauge invariance is preserved. This enables the calcula- 
tion of the effective Lagrangian in terms of the Polyakov 
loop and meson pseudoscalar mesons as basic and inde- 
pendent degrees of freedom. In Sect.[V]we go further and 
integrate over the Polyakov loop variable in a gauge in- 
variant manner, providing the form of the effective (tree 
level) Lagrangian of ChPT at finite temperature. Finally, 
in Scct.|Vl|we summarize our conclusions. 



II. POLYAKOV-CHIRAL QUARK MODELS 

In this section we review the coupling of the Polyakov 
loop to chiral quark models as a quantum and indepen- 
dent variable as suggested in previous work [T^, [20, [2j 
|22|. We follow for the basic ingredients involved in 
the construction of such model. Since the coupling to 
quarks is rather universal we will restrict ourselves to 
the simplest chiral model, namely, a constituent quark 
model (CQ). The corresponding Lagrangian reads 1 



£cq = 9i*9+ i ! + 
=-. qiDq, 



if 



75 - Aft/ 75 - M )q 



(2.1) 

where q = (u,d, s, . . .) represents a quark spinor with 
N c colors and Nf flavours. Mq = diag(rn„, m,i, m s , . . .) 
stands for the current quark mass matrix. The symbols 
(w/,a£) denote external vector and axial-vector fields, 
in flavour space. M is the constituent quark mass and 
U = e 4 ^*//" (y w being the pion weak decay constant 
in the chiral limit) is the flavour matrix representing the 
pseudoscalar octet of mesons in the nonlinear represen- 
tation: 



$ = 







7T+ 


K+ 


7T 








K- 




K° 


-A* 



(2.2) 



For vanishing current quark masses £cq is invariant un- 
der local U (N f) n<g)U (N ~f) l transformations. In addition 
there is a global SU(iV c ) symmetry. 

After formally integrating out the quarks one gets the 
effective action 



r CQ = -zTrlog(*D) 



(2.3) 



1 In Minkowski space we will use Bjorken-Drell conventions 
throughout the paper. 
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We use Tr for the full functional trace, tit stands for 
the trace in flavour space, and tr c for the trace in colour 
space. The ultraviolate divergences introduced by the 
functional determinant can be conveniently handled by 
the Pauli-Villars method. Let us note that the issue of 
regularization is not so crucial in the present case, T <C A 
|l6j . since the divergences affect only the zero tempera- 
ture contributions |33j, |3J] . 

The chiral quark model coupled to the Polyakov loop 
corresponds to introducing a non trivial colour compo- 
nent of the vector field in the temporal direction 



(2.4) 



in the Dirac operator, eq. Q2.ip . The field Vq(x) acts as a 
chemical potential in color space |22|, Upon Wick rota- 
tion to pass to imaginary time (which we use to introduce 
finite temperature) this field gives rise to the Polyakov 
loop 



fi(x, £4) = Texp i 



l-Xi+0 



dx 4 v 4 {x } £4) 



(2.5) 



As discussed in the integration over the Polyakov 
loop coupled to quarks suppresses unwanted coloured 
quark states and mimics the confinement mechanism 
within the quark model. Note that the Polyakov loop 
measure Dfl e ira ^ preserves center symmetry, i.e., in- 
variance under 't Hooft transformations [3^| under which 
fl — > zfi, with z Nc = 1. This symmetry is preserved by 
hadrons but not by quarks. In the Polyakov-chiral quark 
model picture, as a valence quark propagates at finite 
temperature, one can at any place insert a path starting 
and ending at same point x which winds a number n of 
times around the thermal cylinder, picking up a factor 
((— f2) n ) that goes with the propagator factor e~'"' M / T 
(modulo we subleading polynomial corrections). In the 
presence of sea quarks there are further thermal wind- 
ings, giving roughly (being L the number of sea quark 
loops) 



n+ni+- 



■+n L \ e -(\n\+\m\+-+\nL\)M/T 



(2.8) 



where T indicates the Euclidean temporal ordering and 
f3 = 1/T, The Polyakov loop plays an important role as 
an order parameter for the deconfincmcnt transition. In 
the present context it appears quite naturally, since as 
shown in J2^> |3(| a generic effective action with gauge 
fields at finite temperature will depend not only on the 
standard zero temperature local gauge covariant opera- 
tors O n (x) constructed with the covariant derivative, but 
also on the Polyakov loop Cl(x) as a new finite temper- 
ature gauge covariant operator. Specifically, at the one 
loop level the effective Lagrangian takes the generic form 

m 



£{x) =5^tr[/ n (n)O n ]. 



(2.6) 



In the Polyakov-chiral quark model the partition func- 
tion takes the form 



(2.7) 



where DU is the Haar measure over the flavour group 
SV(Nf) (actually, the product of local Haar measures) 
and -Dfi is essentially the Haar measure of the colour 
group SU(iV c ) (see however Sect. 01, Tq is the effective 
gluon action whereas Tcq stands for the quark effective 
action in l|2.3|) . The model is motivated by the underlying 
idea that all zero temperature gluon degrees of freedom 
have been integrated out to yield the constituent quark 
mass, leaving unintegrated the Polyakov loop as the only 
specifically finite temperature gluonic degree of freedom. 
Ideally Tg^] would be the result of such a partial inte- 
gration in pure gluodynamics. The Polyakov loop inte- 
gration implements the gauge invariant integration over 
the w§ field |37|. Since gluons are vector fields, no colour 
axial- vector component is introduced in (|2.4|) . 



The average over the gauge group is that for a pure 
Yang-Mills theory and so preserves triality in the con- 
fined phase that we are considering, i.e. , n + n\ + • • ■ + til 
must be a multiple of N c . For N c > 3, the leading ther- 
mal corrections comes from contributions of the form 
n = —n\ = I, ri2 = . . . = riL = 0. This describes a 
thermal qq pair, with a propag ator factor er 2M l T . Tak- 
ing into account quark interactions the leading thermal 
corrections will be of the type e~ mir l T , being the pion 
the lightest meson. 

A remaining technical point deserves comment. As 
noted above the presence of the Polyakov loop in the ef- 
fective action is an inescapable consequence of having a 
gauge theory at finite temperature. So in addition to 
the standard colour Polyakov loop, there will be also 
Polyakov loops associated to the vector and axial-vector 
flavour external fields. Such chiral Polyakov loops take 
the form 



Qr,l(x) = Texp I i I 1 



dx\ (v{(x, x' 4 ) ±a{{x, x' 4 ))j , 
(2.9) 

and they appear automatically in any gauge invariant one 
loop computation at finite temperature. Such operators 
(as well as the vector flavour Polyakov loop) would in- 
troduce interesting new effects such as center symmetry 
in flavour space, preserved by mesons but not necessarily 
by baryons, depending on the number of flavours. Since 
spin 1 mesons are relatively heavy we do not expect the 
flavour Polyakov loops to be essential in the low temper- 
ature regime to be study in this work, so it will be dis- 
regarded in what follows. Nevertheless they could play a 
more active role near the transition temperature. 
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III. THE STRUCTURE OF THE LOW ENERGY 
CHIRAL LAGRANGIAN AT FINITE 
TEMPERATURE 

As already shown in previous works at zero temper- 
ature |M EE El E3 ( see e -S- for an updated list 
of references) chiral quark models may provide a quan- 
titative and microscopic understanding of the structure 
of the low energy effective Lagrangian of ChPT for the 
pseudoscalar mesons at tree level, namely providing nu- 
merical values for the leading N c contributions to the low 
energy constants (LEC's). In such a framework, external 
currents are minimally coupled at the level of the more 
microscopic quark degrees of freedom. On top of this, me- 
son loops would provide subleading 1 /N c contributions to 
the LEC's in addition to the standard unitarity correc- 
tions of ChPT H3. 

In this section we extend the zero temperature results 
to the finite temperature case and consider also the influ- 
ence of the Polyakov loop in the scheme described in the 
previous section. More explicitly, the partition function 
in 1)2.7(1 can be written as 

Z = J DUDne iV G [U] e i! d^xCnix) ( 31 ) 

and then 

Z = J DU d*xC(x) _ ( 3 2 ) 

The Lagrangian Cq(x) is obtained after integration of 
the quarks, and depends on the fields U(x) and Vi{x). 
Setting Q, = 1 one obtains the ChPT Lagrangian at fi- 
nite temperature corresponding to the CQ model without 
coupling to the Polyakov loop. In the general case Cn(x) 
provides also the coupling of Goldstone mesons to the 
Polyakov loop field, and interesting subject by itself. It 
is not clear, however, how such a Lagrangian could be 
formulated in a model independent way since in QCD 
the Polyakov loop is renormalized and in general will lie 
outside the SU(iV c ) manifold. 



The computation of Cn(x) can be carried out by fol- 
lowing the methods developed in Ref. [3^, |3|| and al- 
ready applied to QCD in [45| • The procedure is detailed 
in the next Section. It is worth noticing that Cq(x), 
either directly, after setting fi = 1, or averaged over 
the Polyakov loop, represents a one quark-loop result 
and hence a quenched approximation to the chiral La- 
grangian. 

Integration over the Polyakov loop field in l|3.1|) gives 
the (unqucnchcd) finite temperature ChPT Lagrangian 
of the model, denoted C*(x), which is the main purpose 
of this work. In ChPT at finite temperature it is usu- 
ally assumed that the corresponding low energy constants 
are temperature independent [f| (see also Refs. [HI. l46j ). 
This is a quite natural assumption because the appli- 
cability of ChPT is based on the existence of a mass 
gap between the Goldstone bosons and the rest of the 
hadronic spectrum. For non-strange mesons the gap is 
provided by the p meson mass My, so one expects the 
temperature dependence of the LEC's to be of the order 
of e~ Mv l T . In a chiral quark model, however, the pseu- 
doscalar mesons are composite particles of constituent 
quarks of a mass M, and the finite temperature also in- 
fluences their microscopic quark substructure. As a con- 
sequence the LEC's become temperature dependent of 
the order e~ 2M ^ T . Our calculation below makes these 
remarks quantitative and also provides an understand- 
ing on how the Polyakov loop cooling mechanism works 
in favor of the ChPT expectations at finite temperature. 

Although the calculation is straightforward, reducing 
effectively to evaluation of Dirac and flavour traces, it is 
technically involved so we discuss here the structure of 
the final result for the general case. More details will be 
elaborated in the next sections. 

The low energy effective Lagrangian written in the 
Gasser-Leutwyler notation and in Minkowski space 
reads as follows (we often use an asterisk as upperscript 
for finite temperature quantities), 



C* {Q) (x) = B* , (3.3) 

f*2 

£*(2)( x ) = ±2_ tr/ (B^WU + x^U + xU^) , (3.4) 

= Ll(tT f (B^U^U)) 2 + L* 2 (tr f (D^r>„U)) 2 

+L;tr / (D M J7 t D^C/D t/ C/ t D ,y C/) + L3tr / (Dof/ t DoC/D AI t/ t D A1 [/) 
+ZJtr / (D M CTtD"l7)tr / 0c t l7 + x rf) 

+L* 5 ti- f CD^T>»U( X ^U + lf x )) +i;tr / (D C/ t Do(7(x t f/ + U* X )) +t*tT f (xD T> U^ + X ^o^ U) 
+L*(tr J ( X t[/ + x^)) 2 + L*(tx f &U - X U r )) 2 
+L'*tr / (f/ t D D C/ - UB -D U^tr f (x^U - X U f ) 
+L*tr / (x t f/x t E/ + X^W) 



5 



-iLlixf{F^Di*U^T> v U + F^WUWU^) 

-ilgtr/ (EfiparfTtU - T> l U^T) U) + £f (DoUTW* - WUT)^)) 

+L* w t If (U^F^UF^ R ) + H*tv f ((F^) 2 + (FL v f)+Hltr f ((E*f + (E^) 2 ) + H;tr f tfx) ■ (3-5) 

I 



Here, tif means flavour trace and we have introduced 
the standard chiral covariant derivatives and gauge field 
strength tensors, 

f; v = i[D*, di] = d,v; - d v v; - i[v^, f;], (3.6) 

with r = L,R, and V^ ,L = ± a?. Finally, as at 
zero temperature, we have normalized the explicit chiral 
breaking field \ = 2i?gM so that £*( 2 ) takes a standard 
form. This introduces a LEC B* {) = \(qq)*\/f* 2 , where 
(qq)* is the quark condensate for one flavour at finite 
temperature. 

The general form of the chiral low energy effective 
Lagrangian, Eq. I|3.5fl requires some remarks. As we 
see there are terms which can be written as the zero 
temperature Lagrangian but with temperature depen- 
dent effective couplings. In addition, we also have new 
terms, which due to the finite temperature generated by 
a heat bath at rest necessarily break Lorentz invariance. 
This feature was already pointed out in 36]. The re- 
markable, not yet understood, feature is that there ap- 
pear less terms breaking the Lorentz symmetry than one 
might naively suggest (for instance a term of the type 
tr/(D £^Dof7) is missing). Thus, some accidental sym- 
metry may be at work and it would be interesting to find 
it explicitly. 

Although the form of the Lagrangian is quite gene- 
ral, the particular values of the low energy coefficients 
obviously depend on the specific model. We will consider 
two particular chiral quark models: the constituent quark 
model and the spectral quark model ^3] • The results for 
the latter model will be presented in appendix [51 and we 
concentrate on the CQ model in the following sections. 

IV. HEAT KERNEL EXPANSION AT FINITE 
TEMPERATURE IN THE PRESENCE OF 
POLYAKOV LOOPS 

The calculation of the effective chiral Lagrangian of 
Eqs. H3.3fl - H3.5fl can be divided into several steps. First, 
we construct a Klein-Gordon operator out of the Dirac 
operator and its adjoint for the non-anomalous part of 
the effective action. After a generalized proper time re- 
presentation we are naturally lead to a heat kernel of 
the Klein-Gordon operator, for which a heat kernel ex- 
pansion is particularly suited. This identification allows 
to directly apply the results of [3(| |45| and to work out 
the traces after applying some relevant matrix identities 



which reduce the number of independent operators. Fi- 
nally, the equations of motion are also considered to take 
into account that the pion fields are on the mass shell. 
In the present section we assume a fixed value of the 
Polyakov loop. This provides £n(x) of 1)3. 1|) . In Sect. E] 
the Polyakov loop integration will be addressed, to yield 
the unquenched chiral Lagrangian C* (x) of i|3.2|l . 

In this section and the next one we revert to Euclidean 
space, X4 = ixo, since it is much more convenient for 
calculations at finite temperature in the imaginary time 
formalism. The partition function is numerically un- 
changed, but takes the form 

Z = J DUe-$ dlxC ' {x) (Euclidean). (4.1) 

Likewise, the low energy coefficients are also unchanged. 
E.g., 

C* (tt) {x) = -B* (Euclidean) (4.2) 
C* {2) {x) = £.tr f (DrrtT> li U-x*U- X rt) , 
with the same values of B* , /*, etc. 



A. The effective Klein-Gordon operator 

The Dirac operator appearing in the fermion determi- 
nant behaves covariantly under chiral transformations. 
This implies that, in principle, one has to consider both 
vector and axial- vector couplings. A great deal of sim- 
plification is achieved if the conventions of Ref. [i^ . |49| 
are considered, where it is shown that it suffices to carry 
out the calculation in the simpler case of a vector-like 
coupling and then reconstruct the total chiral invariant 
result in a suitable way. 

Let us hence consider the following Dirac operator with 
a vector-like coupling, 2 

B=p+h, h = m + z, (4.3) 

where h includes the pion field to, which we take 0(p ) 
in the chiral counting, and the chiral symmetry breaking 
mass term z which we take 0(p 2 ). We also introduce the 



2 Dirac gammas in Euclidean space are taken Hermitian, and 
7n7^ = S^iu + o" M „. 
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following useful notation 

m LR = MU , m RL = MC/t , (4.4) 

1 1 + , A . 

zlr. = i^X, z RL = -—x\ (4.5) 

where \ has been introduced after l|3.6l) . (In the notation 
of |48j,|49j, the symbol m is to be interpreted as tulr or 
mRL depending on its position in the formula, and sim- 
ilarly for the other vector- like symbols.) It is convenient 
to introduce the adjoint Dirac operator 

D f = -#>+/!,. (4.6) 

This definition allows us to separate the effective action 
into 75-even and 75-odd components, corresponding to 
normal and abnormal parity processes, since they corre- 
spond to the real and imaginary parts of the Euclidean 
effective action, respectively (see e.g. |50j). 

We will focus on the normal parity component of the 
effective action, which from Eq. (|2.3() is formally given 
by 

T+ Q = -iTrlog(DtD) =: J dx 4 J d 3 X £ n (x) , (4.7) 

with the relevant Klein-Gordon operator 

D+D = -Dl + ^a^F^ - 7 „D„/i + z 2 + {m, z] + M 2 , 

(4.8) 



We use the notation D M /i = [D^,h]. The field strength 
tensor is defined as F^ = i[D^, D v ]. The operator D^D 
is of the Klein-Gordon type, with mass term 

^(J^uFy.u - 7^/1 + 7? + {to, z} + M 2 . (4.9) 

Therefore a heat kernel expansion becomes appropriate. 
After a (generalized) proper-time regularization, the ef- 
fective Lagrangian in Euclidean space becomes 

C n {x) = - / — 0rtr(x|e- TDtD |x) (4.10) 
2 Jo T 

1 r°° dr p-rM 2 

2 Jo T (4ttt) 2 ^ 

This proper-time representation has been chosen to ac- 
commodate the Pauli-Villars regularization as used in 
e.g. Ref. 0, 

<H r ) = J2 c * e " TA ' > co = l,Ao = 0, (4.11) 

i 

but still using full advantage of the heat kernel. The 
conditions J2i c i^i = 0; f° r n ~ 0,2,4, allow to render 
finite the logarithmic, quadratic and quartic divergencies, 
respectively. To 0(p A ) terms we obtain the following 
contributions for the thermal Seeley-DeWitt coefficients 
36, 45] after the Dirac trace has been implemented, 



u l/2 

K 

t>3/2 

K 

K/2 

h* 
"3 

t>7/2 

b*A 



4^0 (n) , 
0, 

-4^o(^)({to,z} + z 2 ) , 
0, 

2po(fi) ((m M ) 2 + {m^ z^} + {to, z} 2 + ±F 2 ^j + j^ a (fi)£? + 0(p 6 ) , 

o( P 5 ), 

2 



(fio(tl) ( m M {m M , {to, z}} + {to, z}(m^) 2 - £{-F/«» "fym„} + im ll F ilu m v + -(m M „) 2 



(4.12) 



In these formulas Ei — Fa is the "electric field" and 
a notation of the type X^ va has been used to mean 
D M D„D Q A = [Dfi, [D v , [D a ,X]]}, e.g. to 4m = D 4 D m to, 
Fafj, v — D Q F M „. The functions ipo and Tp 2 are defined in 
appendix lAl 

Note that, in principle, since -D M contains both flavour 



and colour gauge fields, the Polyakov loop and the field 
strength tensor F^ u (and derivatives) in the heat kernel 
coefficients will also contain flavour and colour contribu- 
tions. As explained before in the Polyakov loop we retain 
only its colour part. In addition, in the strength tensor 
and derivatives we keep only the flavour part. This ne- 
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gleets some gluonic corrections. This kind of QCD cor- 
rections have been considered for instance in |4lj for the 
Nambu-Jona-Lasinio (NJL) model. We disregard them 
in the present treatment since they are not specific of 
finite temperature. 



B. Effective Lagrangian 

The effective Lagrangian can be written as 

Cn(x) = £$\x) + Cg\x) + C^(x) +-.., (4.13) 

where Cffl is of 0{p n ). Making use of the expression 
(|4.10() , the flavour trace of the Seeley-DeWitt coefficients 
obtained in (|A10f) and after evaluating the proper-time 
integrals using the Pauli-Villars regularization (see ap- 
pendix^), we get for the Lagrangian Cq(x) an expres- 
sion formally identical to (the Euclidean version of) that 
in (|3 . 313. 5|) . except that the coefficients depend on the 
Polyakov loop, that is, B*(Q.) instead of B* , etc. 
At zeroth order one finds 



B*(0) 



The second order gives 



2N f M i ^ . . 
-^tr c T_ 4 (fi). 



M 2 



Air 2 



(4.14) 



(4.15) 



Note that we have not yet averaged over the Polyakov 
loop and the (Polyakov loop independent) normalization 
constant B$ (needed to fix the definition of \) is deter- 
mined after this average, as we will explain in the next 
section. However the combination x/2£?o = Mq is un- 
ambiguous. It contributes to (Euclidean) with 



(An) 2 B* 

For the fourth order one obtains 
1 



tr c J_ 2 (r»)tr / ( X t(7 + [/ t x)j (4 . 16) 



L$(n) 

L* 5 (Q) 
L* 6 (Q) 



tr c 2 4 (n) 



24(4tt) 2 
-8Ll(n) + ~L* 9 (n) 



W trcJ2(Q) 



0. 



M //; 2 (0) 



2BX 



AM 2 



3Lj(n) 



L' 5 *(Q) = 



= o. 



8N f 



n 2 m 

2B*M 



4.17) 
4.18) 
4.19) 

4.20) 
4.21) 
4.22) 

4.23) 
4.24) 

4.25) 



L (ft) 
Lt(Sl) 

L* 9 (n) 

Km 



AN f 
1 ( 1 



165* \M B* 



L' 9 *(n) = -T 3 (ri), 
i 

r 



■o £ 9( fi )> 



L* g (n). 



2AM 2 A 

8Br 



1 



2 ■ 



(4.26) 

f: 2 (n) - ^L5(nj4-27) 

(4.28) 

(4.29) 
(4.30) 

(4.31) 
(4.32) 

(4.33) 



Note that all new Lorentz breaking terms, except H 1 , 
are proportional among them. On the other hand, the 
standard Gasser-Leutwyler coefficients can be expressed 
in terms of /* 2 ,i?o,L^ and L 9 or, equivalently, in terms 
of the integrals ti c I n for n = —2, 0, 2, 4, which are com- 
puted in appendix 1X1 



V. POLYAKOV LOOP INTEGRATION 

A. Colour Group averaging 

In order to proceed to the computation of the chi- 
ral (Euclidean) Lagrangian we need to carry out the 
Polyakov loop integration, 



fd 4 xC'(x) = / Dne -Ta[n] e - J d 4 xCu(x) 



(5.1) 



To this end we can take advantage of the chiral counting 
to expand the exponential with Cq(x) to fourth order. 
There is an obstruction, however, since Cq^ will appear 
to all orders, making it difficult to take the Polyakov 
loop average. To sort this problem we introduce a fur- 
ther thermal counting, in addition to the chiral one. The 
thermal counting suppresses terms which become irrel- 
evant at low temperatures and so it is compatible with 
the chiral expansion. 

The thermal counting is defined as follows. Cn(x) 
comes from one quark-loop. As noted at the end of 
section [H] any such loop will wind £ times around the 
thermal cylinder (in imaginary time). Paths with I = 
give the zero temperature Lagrangian Co which depends 
on U but not on f2. The other paths pick up a factor 
tr c [(— fi) ] which is accompanied by a suppression factor 
(at low temperature) of the order of e~\ t -\ M l T from the 
quark propagator. They give the thermal component of 
the Lagrangian. 



Cq(x) = C (x) + £ei x ) ■ 



(5.2) 



8 



In this thermal counting the terms are ordered by the 
number I of Polyakov loop they carry with tr c [(— S7) £ ] as 
order \t\. The discussion of the incidence of the gluonic 
part of the action Yq [O] on this thermal expansion will 
be postponed till Sec. IV Dl 

Combining the chiral and thermal expansions, one has 



(5.3) 



n.f 



(Note that the chiral expansion of the zero temperature 
component Cq(x) is not required since it does not depend 
on the Polyakov loop.) This is to be introduced in the 
Boltzmann factor exponential of l|5.1(l and expanded. We 
will retain terms of 0(p 4 ) in the chiral expansion and of 
leading order in the thermal expansion. Because the Haar 
measure Dfl (and actually also the gluonic correction) 
preserves center symmetry, the first thermal correction 
with a single quark loop (i.e., keeping just one Lagrangian 
in the expansion of the exponential) will be of order £ = 
N c (a baryon-loop like contribution). For N c > 3 this 
is not the dominant contribution. The leading thermal 
correction comes instead from correlation of a quark loop 
with an anti quark loop, a meson-loop like configuration. 
Schematically, these corrections have thus the following 
structure 



£(0)^(0) + £(°)£( 2 ) 



4 0) 4 4) 



£(2) £(2) 



(5.4) 



contributing to £*(°), C*^ and respectively. Note 

that the two quark loops occur at different spatial points, 
so the correlation between two Polyakov loops is needed. 
This we model as Il3 



Df]tr c fi(f)tr c fi _1 (?7) = e 



-a\x-y\/T 



(5.5) 



where a = (425 MeV) 2 is the string tension. (This, of 
course, implies that Dil is not exactly the product of local 
Haar measures.) This prescription is consistent with the 
zero temperature quark-anti quark potential at medium 
and large distances, 3 while at coincident points it repro- 
duces the SU(iV c ) group integration formula |5l). 



dQ tr-fi tr„f2 



1 . 



(5.6) 



Effectively, the consequence of using (|5.5(1 in a two-point 
computation of type l|5.4|l is to introduce a correlation 
domain of spatial size V 

J d 4 xd 4 y J DntT c n(x)tr c n-\y) = J d A x^ , 
8ttT 3 



V = 



(5.7) 



3 The short distance Coulombian part of the qq potential cannot 
be reproduced with unrenormalized Polyakov loops lying on the 
group manifold. 



Application of the procedure just described produces, 
for the zeroth order of the effective Lagrangian in the 
chiral expansion, 



B* = B 



NfM 4 T 3 V 



-Ki 



(5.8) 



In this formula K 2 = K^jM/T) and K n (z) refers to the 
n-th order Bessel function [52]. We use the same con- 
vention in all forthcoming expressions for the LEC. Not 
surprisingly, the thermal correction to the vacuum en- 
ergy density is proportional to the correlation volume V 
and contains exactly two Bessel functions from the two 
correlated quark loops. This holds also for all LEC to 
follow. 

To 0(p 2 ) we find 



f* 2 = f 

J IT J TV 



B* 



Bo 



N f M 4 TV 
NfM^TV 



K Q K 2 , (5.9) 



K 2 {B0K0 - 2TK X ) .(5.10) 



And therefore, for the single flavor quark condensate at 
finite temperature 

2N f M i T 2 V 
(qq)* = (W + f 4 K1K2 (5.11) 

At very low temperatures, the asymptotic form of the 
Bessel function 



sempe 

m 

*«w - (£) 



1/2 



(5.12) 



applies, and so the first thermal correction in the low en- 
ergy coefficients has the exponential suppression e~ 2M l T . 
This indicates that meson thermal loops dominate at low 
temperatures (in our treatment we ignore completely the 
quark binding effects so the nature of the mesons involved 
cannot be resolved within the present approach). Note 
that the group integration produces an A^ c suppression 
with respect to the zero temperature contribution [Tol]. 4 
We warn, however, that the results obtained with the full 
Bessel functions and with their asymptotic form differ be- 
yond very low temperatures, especially as n increases. 

Due to center symmetry of the Polyakov loop measure, 
the higher thermal contributions are of the general type 
J DfltT c (n ni ) ■ ■ -tr c (fT f ), with a multiple of N c . 

Among other, there are 

• Two or more meson-like loops, of the type 
/Df^trcfHrc^- 1 )", n = 2,3, . . ., with a thermal 
suppression factors e - 2nM / T . 



4 Naively, a two quark-loop term would be O(N^), and in fact, 
ignoring the Polyakov loop dependence, i.e., setting Q = 1 as 
in the standard NJL model, would give , instead of unity 
on the r.h.s. of 15.61 . Thus the introduction of the Polyakov 
loop proves essential to reproduce the 0(N®) result that one 
should expect for the thermal corrections coming from colourless 
hadronic excitations. 
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• Baryon-like loop excitations, of the type 
J DQtr c Q ±Nc , with thermal suppression e~ NcM ' T '. 
They come from a single quark looping iV c times 
the thermal cylinder. These are the leading ther- 
mal corrections of hadronic type in the quenched 
version of the theory [53l . 

• Gluonic corrections, with thermal suppres- 



sion e 



-ma/T 



, m G ~ 650 MeV, from T G [Q] 



In passing, we note that the leading thermal contri- 
bution to the Polyakov loop expectation value can also 
be easily computed. It comes from the coupling of the 
Polyakov loop observable with a quark loop. Schemati- 
cally, from (flC^). An easy computation gives, at lead- 
ing order 



,1 n , N f M 2 TV „ 



(5.13) 



This is a measure of the explicit breaking of center sym- 
metry by quarks, and is controlled by the ratio M/T in 
the exponential factor. 

B. Low Energy Coefficients 

Let us come now to 0(p 4 ) terms. The result we find 
for the Gasser-Leutwyler coefficients for the constituent 
quark model reads 



T* T + M * V (k 2 N fK* 



NfM 4 V 2 



(5.14) 
(5.15) 



2 192tt 4 T 
NtM 3 V 

L l = L 3 + / 87r4r K * ( MK * - TK ^) > ( 5 - 16 ) 

— * N f M 4 V 

3 487r 4 T ° 2 ' v ; 

(5.18) 



M ' V K K 



La = T~ 

NtM 3 V 

L* 5 = L 5 - eL 5 + ^ K 2 (MKi - 2TK ^.19) 

L* 5 = t* = , (5.20) 
M 4 TV , 

L « = 64^ < (5 ' 21) 
f 2 

M*V /12T \ 

192tt 4 V B o ' 



L = - 

47V. 



(5.23) 



8 I6B0 V^o M/ 192tt 4 



xlKl + *«L(±-±i K . ] , (5 . 24) 



B V^o M 



^9 - ^9 



N f M 3 V 
24tt 4 
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= # 1 



N f M 2 V 



24tt 4 



^1^2 , (5.25) 

(5.26) 
(5.27) 

if 2 (4TiY - MKi) , (5.28) 
(5.29) 



Hr, = H2 



f 2 

N f M 3 V „ (Y1MT 



-K- 



96^ 4 4 \ B 2 
where we have defined the quantity 
e=(B*- B )/Bo . 



Ko-KA ,(5.30) 



(5.31) 



As usual, to bring the fourth order Lagrangian to the 
standard form 1|3.5[) we have used the equation of mo- 
tion for the nonlinear U field. This is obtained by min- 
imizing C*( 2 \x) (and hence after the average over the 
Polyakov loop) and reads 



B* [u,z} + ^tr f ([u,z]) = 0. (5.32) 



(where the vector-like notation has been used, with 
ulr — U and url = U*). The last term arises from 
the constraint det(C7) = 1, since we are considering a 
SV(Nf) flavour group. 

Explicit expressions for the zero temperature coeffi- 
cients appear in Once again, we observe that the 
finite temperature corrections are iV c -suppressed as com- 
pared to the zero temperature values, as expected from 

hadronic excitations, and this requires a proper Polyakov 

— '* 

loop average. The terms with Nt, and also L£ and L , 
are those coming from while those without come 

from Ci Ci . It is worth noticing that some struc- 
tures which were absent in the quenched approximation 
Cq\x), e.g. L\ and Lg, are allowed in the unquenched 
result, from C± £± . 



C. Volume rule 

Clearly, to carry out the functional integration over 
the Polyakov loop indicated in i|5.1[l , besides requiring a 
careful specification of the Polyakov loop action, would 
be an exceedingly demanding task. However, we can 
make systematic relations of the type l|5.7fl using a simple 
model. This will allow us to go beyond the leading ther- 
mal corrections and to include gluonic corrections from 
Tq [Q] . Specifically, we assume that the space-time is de- 
composed into correlation domains of size V/T, in such 
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a way that two Polyakov loops are completely correlated 
if they lie within the same domain and are completely 
decorrelated otherwise. So for a partition function of the 
form 

Z = J Dne-$ d * xC ^ :=e ~.fd 4 xC( x ) ( 5 33 ) 
we take 

Z= I Y[dn n e-^^ r " (5.34) 

J n 

where n labels the correlation domain, dfl is the Haar 
measure, and T n = (V/T)C(x n ), x n lying in the domain 
ii. This gives 

z= [Y[dn n j i-^r n + ~^r„rv + ••• 

n \ n n,n' J 

= i - £<r n ) + \ j^K) + < r »>< r «'> + ■■■ 

n n n^n' 

(5.35) 

where ( ) indicates average of fl over SU(-/V C ). Obvious 
rearrangement gives, finally 

£{x) = (£(*)) - ~ ((C 2 (x)) - (C{x)f) + ■■■ (5.36) 

This is the standard cumulant expansion. Equivalently, 
C'(x) can be read off from 

e -WW = ( e -WA (5 37) 

\ /SU(AT C ) 

In concrete configurations the r.h.s. can be computed 
numerically (being the integration manifold very well be- 
haved.) Points to be noted are: i) The effective action 
is an extensive quantity. It is given by J d A x£(x) and 
not (y/T)£ (x). ii) If C{x) contains a piece C {x) which 
is Polyakov loop independent, this term simply adds in 
C'(x). It does not appear in the ^-dependent contribu- 
tions. And, iii) previous results of this section, e.g. Ij5.7jl . 
are correctly reproduced. 

D. Gluonic Corrections 

Up to now in l|5.1|l we have neglected the factor e~ r<3 
in the Polyakov loop measure. Following we adopt a 
simple model inspired on the lattice strong coupling ex- 
pansion I55L l56^| . This model has one parameter which is 
fitted in [22J to reproduce the transition temperature in 
gluodynamics in the mean field approximation. (A more 
refined fit involving more terms is discussed in [23.) Be- 
cause our calculation does not rely on the mean field ap- 
proach, our implementation of the model is not identical 
to that in those works and so refitting would probably 
be required. Nevertheless, this will not be needed here 



to make low temperature estimates since the gluonic cor- 
rections are much suppressed. Specifically, since we want 
to have extensivity, we introduce a local Lagrangian 

r G [to] = J d 4 xC g (x) (5.38) 

which is modeled as 

C g (x) = -^Te-^Mft)! 2 (5.39) 

with an adjustable parameter a -1 = 272 MeV. Triality is 
preserved, but, unlike the pure Haar measure, at higher 
temperatures the action favors center of the group val- 
ues for the Polyakov loop. (The spontaneous breaking 
of center symmetry in gluodynamics indicates that a La- 
grangian with a different form would be needed in the 
high temperature regime.) At low temperature we can 
see an exponential thermal suppression controlled by a 
gluonic mass mc — ca = 664 MeV. 

The gluonic corrections can be accounted for by includ- 
ing C g (x) in the Lagrangian of the volume rule l|5.37|l . 
The leading gluonic correction is 0(e~ mG ^ T ) and goes to 
the vacuum energy density. Using l|5.6[) un easily finds 
S g B* = -{C g [x)) = 6a- 3 Te~ mG / T . At low temperature 
this would be comparable with the meson-loop like terms, 
however, recall that the constituent quark mass M (ac- 
tually M*) gets reduced as the temperature approaches 
the chiral transition, thereby enhancing the role played 
by mesonic loops. 

The gluonic corrections to /* 2 , Bq and the other low 
energy coefficients L* are of £)(e~( 2M+mG ^ T ), as is easily 
verified, thus they can be neglected as compared with the 
hadronic-loop corrections. 

Instead of modeling gluonic corrections through a local 
Lagrangian, another approach suggested by the volume 
rule (|5.37|) is to remove the Polyakov loop action from 
the Lagrangian and instead, to replace the Haar mea- 
sure average over the colour group by a weighted aver- 
age, dO, — > dfle~ r °^ . In principle such a weight could 
be obtained through a proper sampling of the Polyakov 
loop distribution on the group manifold within lattice 
computations of gluodynamics. At present this is not 
feasible due to severe numerical and conceptual prob- 
lems related to renormalization issues. An advantage of 
this approach would be that expectation values of the 
type (f(tt(x))} (a single point) would be exact when the 
quarks are switched off. However, it misses the normal- 
ization of the measure, which is needed for its contri- 
bution to the vacuum energy density. From (|5.37|l it is 
clear that such weight is equivalent to a local Lagrangian 
C g {x) = (T/V)T G . 

VI. CONCLUSIONS 

The previous results clarify the question that chi- 
ral quark model practitioners might legitimately ask, 
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namely, whether the tree level low energy coefficients in 
the chiral Lagrangian do genuinely depend on tempera- 
ture or not. According to our present calculation this 
temperature dependence at very small temperatures is 
tiny. It is exponentially suppressed by a scale which is 
exactly the mass gap, and becomes inessential on temper- 
ature scales below the deconfining transition (the region 
of temperatures where the Polyakov cooling mechanism 
acts) if the Polyakov loop is first introduced and colour 
singlet projection is subsequently carried out in a gauge 
invariant manner. In a quark model without Polyakov 
loop the temperature dependence of the low energy con- 
stants behaves as L* — Li ~ N c e~ M / J ', due to spuri- 
ous contributions of colour non singlet multi-quark states 
while the Polyakov cooling, i.e., the explicit suppression 
of colourful excitations provides instead the behaviour 
L* — Li ~ e -2M/T ag ex p ec ted from ChPT arguments. 
Remarkably, this cooling mechanism produces also the 
correct N c counting. As shown in previous work piil | this 
effect at low temperatures manifests quantitatively in the 
chiral symmetry restoration-center symmetry breaking 
phase transition. 

To see how the agreement of Polyakov Chiral Quark 
Models to the theoretical assumptions of ChPT at fi- 
nite temperature well below the phase transition materi- 
alizes in practice we have calculated the chiral effective 
Lagrangian at finite temperature. As a by product the in- 
teraction between Polyakov loop and Goldstone bosons 
may be analyzed in some detail. The resulting chiral 
Lagrangian can be decomposed into a zero temperature 
like looking piece with temperature dependent low energy 
constants and a new Lorentz breaking contribution with 
novel structures generated by the heat bath at rest. We 
remind that these calculations aim at describing the tree 
level effective action of ChPT at finite temperature. In 
any case, the corresponding temperature induced effects 



on the low energy constants at this level of approximation 
exhibit the Polyakov cooling. In other words, below the 
phase transition any temperature dependence on the tree 
level low energy constants may be neglected. This is pre- 
cisely the starting assumption of ChPT theory, a purely 
hadronic theory, and as we have shown it arises naturally 
when the local and quantum nature of the Polyakov loop 
is taken into account. 
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APPENDIX A: FLAVOUR TRACES AND 
PROPER-TIME INTEGRALS 

In this appendix we compute the flavour traces of the 
Seeley-DeWitt coefficients of Eq. (|4.12|) . and the proper 
time integrals which appear in Eq. (|4.10() . in order to 
obtain the effective Lagrangian of sec. II V ij| Euclidean 
signature is used throughout. 



1. Flavour traces and useful identities 

For Nf = 3 flavours we have the SU(3) identity 

tr(ABAB) = -2tr(A 2 B 2 ) + itr(A 2 )tr(B 2 ) + {tv(AB)) 2 

(Al) 

where A and B are 3x3 Hcrmitian traceless matrices. 
From here one gets 



tv f ((m^) 2 ) = -2tr / ((m M ) 4 ) + i(tr / ((m / ,) 2 )) 2 + (tr / (m M m l/ )) 2 , (A2) 
tr f {{m 4 mrf) = -2tr / ((m 4 ) 2 (m M ) 2 ) + ^((to^ 2 )^^) 2 ) + (tr^m^)) 2 . (A3) 
Further useful identities are 

tr/((nv,) 2 ) = tr / ((m WI ) 2 ) + 2itr / (F M „m M m„) - tr f {{mF^) 2 ) + M 2 tr f (Fp , (A4) 
tr/((m 4AI ) 2 ) = tr/ (m 44 m w ) + 2itr f (E i [m 4: , mi\) + 2itr/(£ 4i mTO;) , (A5) 

where we have applied the identity X^ v — X u ^ — i[F^,X], being X any operator. Using the equation of motion, 
Eq. H5.32fl . we obtain the following identities 

tr/(m„* M ) = ^pp tr /(("V) 2 "^) " i4M tr/((m - T)2) + 4^ tr/(x2) + 8M N f B* (tr/(K x])f ' (A6) 

1 1 M 2 1 

tr/((m w ) 2 ) = ^tr / ((m Al ) 4 )--tr / ((mx) 2 ) + — tr f (x 2 ) + — (t r/ ([m, z])) 2 , (A7) 

tr f (m 44 m^) = ^tr / ((m 4 ) 2 (m A1 ) 2 ) - Mtr f (mux) - -^tr f ((m 4 ) 2 mx) + 2Aj L ^ tr / (m 44 m)tr / ([m, x}) , (A8) 
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where the normalized field x — 2B'^z has been introduced Applying Eqs. I|A2|) - (|A8|) the flavour trace of the Seeley- 
and so DeWitt can be worked out, yielding 



xlr = X , X RL = T • (A9) 



tr/65 = 4N m (tt), 



tifbl 



ti f b* 4 



^o(fi) ( ^tr/(ma;) + -^tr^a; 2 ) 



2^ (^)tr/((m p 



R*2 
-°0 



BqM 2 



Lpo{tytrf{{rn ll ) 2 mx) 



1 / 1 



V-Bq m 



ip (fl)tTf((mx) 2 ) 



-i^o(!l)tr/(F F m p m„) + ^ (S])trj(K P ) 2 ) - ^M a po(n)tr/(J^,) - ^MVo(0)tr/(a; 2 ) 

12 1 Af 

+ -VJ (fi)tr / ((ma;) 2 ) - -^(po(^)tr f ([m^) 2 mx) - — — ^ 2 (fi)tr / ((m 4 ) 2 ma;) - —7p 2 (n)tr f (m 44 x) 



3M 



+^^ 2 (^)ti7(^[ m 4,"ii]) + qi^^tr/^immi) - ^F2^ (^) tr /(( m i") 4 ) + ^j ( P2(^) tl 7(( TO 4) 2 (m AI ) 2 ) 



1 



127V/ 
1 



^o(n)(tr/([m,a;])) s 



6MN f 



3Af 2 

^ 2 (^)tr/(m44m)tr/([m, a;])) . 
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\ 2 1 2 
^( fl )( tr /(K) ) - ^<Po(^)(tr/(m AI ™ l ,)) 2 + -po(0)tr / ((m fl ) ) . 



(A10) 



2. Proper-time integrals with Polyakov loop 

The basic proper-time integrals are defined by 

dr 



l 2l (A,M,fl) := M 



j 2i (A,M, n) 



Af 



(T)T l e-™' 



^ 2 (O0412) 



where is a SU(iV c ) matrix in colour space, and we define 
<po and Tp 2 as follows 



„2 a 2 

n 2 e"— (-n) T 



(A14) 



The n = term in the sum corresponds to the zero tem- 
perature contribution, and for such a term the Pauli- 
Villars regularization can be applied. In the remaining 
I! / terms the regularization will be removed, since 
these thermal terms are ultraviolet finite. This approxi- 
mation is justified for temperatures well below the cutoff. 
(A13) The calculation of the integral yields 

I 



M 2 l- 2 {A,M,ri) 
T (A,M,Q) 

i 2i (A,M, n) 



1 CO 

-~ £ c 4 (A 2 + M 2 ) 2 log(A 2 + Af 2 ) + 8(A/T) 2 £ 



K 2 (nM/T) 



^ c 4 (A 2 + M 2 ) log(A 2 + Af 2 ) + 4MT £ K ^JI T ) 4 , o , ■„ , 



n=l 



5>log(A 2 + Af 2 ) + 2 £jr (nM/T) ((-«)" + (-0)- 



n=l 



r(OE' 



Af 2 



A 2 + Af 2 



2 (If) E^^( nM / T )((- Q )" + (-")" 



2r 



I oo 



Re(0 > , 
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t 2J (a,m, n) 



M 
2T 



2 + 1 oo 



n i+1 AViMf/T)((-ft)" + (-n)" 
I 



I e 



(A15) 



n=l 



APPENDIX B: RESULTS FOR THE SPECTRAL 
QUARK MODEL 

In this appendix we present the results of the low en- 
ergy coefficients for the Spectral Quark Model (SQM). In 
the SQM the effective action reads 



SQM 



dujp{u)Tr log (iD) . 



where the Dirac operator is given by 



(Bl) 



(B2) 



and p(uj) is the spectral function of a generalized 
Lehmann representation of the quark propagator with 
u> the spectral mass defined on a suitable contour of the 
complex plane The use of certain spectral 

conditions guarantees finiteness of the action. 

A judicious choice of the spectral function based on 
vector meson dominance generates a quark propagator 
with no-poles (analytic confinement). More details of 
the SQM for at zero and finite temperature relevant for 
the present paper are further developed at Ref . [l!j . The 
partition function for the SQM can be written as 



DUDfi e iV ° [n] e irsQM [u ' n] 



(B3) 



In this model one has to compute an average over the 
constituent quark mass with an spectral weight function 
p(u) satisfying a set of spectral conditions Note that M 
appears as argument in the integrals T 2 \ and Z 2 z, but 
also as additional multiplicative factors in the Gasser- 
Leutwyler coefficients. This generates a larger number of 
independent functions as compared to the NJL model. 

We write first the results for the effective Lagrangian 
before the group average. We get for the zero order La- 
grangian 



,(o) 



2Nf 
(4^ 



(w 4 tr r Z_ 



(B4) 



To simplify the notation we indicate with ( ) the spec- 
tral average J dujp(cu). The second order Lagrangian be- 
comes 



(2) _ ff{n) 



(47T)2B * 



tr/(K) 2 ) 

(w 3 tr c J_2(^))ti7(ua;), (B5) 



with 



f; 2 (n) = _<c 2 tr c j (tt)) 



(B6) 



As in the CQ model, Bq will be obtained after the average 
in Q. The fourth order Lagrangian has the low energy 
coefficients 



Lt(Q) = 



1 



L* 5 (Q) = 

L* 7 (n) = 

L* 8 (S1) = 

L* 9 (n) = 



24(4tt) 2 
1 

6(4tt) 
1 



(tr c I 4 (tt)> 



2(4tt) 2 BS 
1 



(cotr c I (n)) - (wtr c :i 2 (fi)) 



1 



2{An) 2 N f V 2B 

(ajtr c 2o(^)} 



4(4^) 2 BS 



16^ 2 8 



L* 9 (Q), 



1 



3(4^) 2 
1 



(tr c x 2 (r>)} 



1 



6(4tt 
1 



o<tr c 2To(0)) 



2(47r) 2 BS \B* 

fz 2 m , i 



— (w 2 tr c 2_ 2 (ft)) - (wtr c J (Q)) 



, , -iS(fi). (B7) 

8BX 2 4 9V ' K ' 



The remaining coefficients satisfy the same geometric re- 
lations obtained in the NJL, Eqs. (|4. LTf) - (|4.33|) . We also 
get the relation 



L* 7 (n) 



N f V 16B * : 



(B8) 



both in the SQM and in NJL. 

The next step is to carry out the integration over 
Polyakov loops, as in Sec. For the temperature in- 
dependent part the result can be expressed in terms of 
the spectral log- moments as was already done in Ref. ^2 . 
For the temperature dependent contributions we will ex- 
plicitly use the meson dominant form of the spectral func- 
tion p(ui). After computing all averages (spectral and 
colour group average), we get for the zeroth order effec- 
tive Lagrangian at leading order in the thermal expansion 



B* = B + 



NjT 3 V 



V^2> 5 



(B9) 



In this expression K 2 = K 2 {uo/T). We use the same 
convention for the symbol K n in subsequent expressions. 
The spectral average can be computed and gives 

T 2 

(uj 2 K 2 ) = —(48 + 24x v + 6a; 2 , + x\)e~ Xv I 2 . (BIO) 
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We use the notation xy — My/T. 

The pion weak decay constant and the normalization 
constant read respectively 



f*2 

J 7T 
ft 

Bo 



1 - 



= 1 



N f TV 



{lj 2 K q }(uj 2 K 2 }, 



^ft 



9T 

x ( (u?K Q )-— (w 2 K x ) 

Bo 



with the spectral averages 
T 2 



(<SKo) 



24 



x 2 v (2 + x v )e~ xv/2 



(Bll) 



(B12) 



(B13) 



(oj 2 K x ) 



£3_ -xs/2 

AT 



(B14) 



and xs = Ms/T. The equations of motion are those in 
Eq. (|5.32|> . (Note that Bq, and hence the normalization 
of the field \ depends on the model.) The low energy 
coefficients, including the first thermal correction, read 



L\ 
LI 
L* 3 

T * 
T * 

T 

L* 
T 



T * 



HI 

Hi 
H* 



Li 
L 2 



2V 



3(4tt) 4 T 

N f V 
3(87r 2 ) 2 T 



[co 2 K ) 2 -N f {co 2 K 2 ) 2 ), 



= u 



NfV ;(oj 2 K 2 )( (lo 2 K 2 ) -T 2 (ojK x ) 



3(2tt) 4 T 
NfV -{co 2 K )(co 2 K 2 ) 



3(2tt) 4 T 



V 



(2tt) 4 B 



(uj 2 Ko)(uj 2 K x ) 



(1 - e)L 6 - 1 ^_(uj 2 K 2 ) (2T(ujKo) - (uK x ) 



L 5 - ^ L 3 , 

TV 
64n 4 B 2 
L 



V 



/12T 



= L* 



8B 2 2 



12T ( 1 

-B \B Q 



L Q = 



L 9 - ^^(uK x )(u; 2 K 2 ) 
t* - 7* 



-L. 



9 • 



= H x 



NfV 
6(2tt) 



{luK x ){lo 2 K 2 ) 



NfV 
24tt 4 



(B15) 
(B16) 
(B17) 
(B18) 
(B19) 
(B20) 
(B21) 
(B22) 

(B23) 

(B24) 

(B25) 

(B26) 
(B27) 
(B28) 

(B29) 
(B30) 



1 9T f 1 

— ((u,K ) + — (u; 2 Ko) 
t>o V y o 
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24T 2 

(u,*Ki)-— 2-(uKi 



(B31) 



where e = (Bq — Bq)/Bq. The spectral averages are 
{Ko) ■ 



.I 7E -log(^/4) + ^(5/2) 



7200 



L F 2 [5/2; 7/2, 7/2; (xy/4) 2 



^ 2 F 3 [1, 1; -1/2, 2, 2; (zy/4) 2 ] (B32) 



P3 

2T 2 x% 



(2 + xs) e ^ s/2 , 



(B33) 



(uK t ) = -f-{12 + Qxv+x 2 v )e- Xv l 2 
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(B34) 



pF q [ai, . . . , ct p ; &i, . . . , b q ; z] is the generalized hypergeo- 
metric function. The expressions for the zero tempera- 
ture coefficients appear in |4^ . 
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